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It has been a big challenge for lattice QCD to simulate dynamical quarks near 
the chiral limit. Theoretically, it is well-known that the naive chiral symmetry 
cannot be realized on the lattice (the Nielsen-Ninomiya theorem). Also practi- 
cally, the computational cost rapidly grows as the quark mass is reduced. The 
JLQCD collaboration started a project to perform simulations with exact but 
modified chiral symmetry using the the overlap-Dirac operator and the topol- 
ogy conserving action. The latter is helpful to reduce the numerical cost of the 
dynamical quarks. Our simulation of two-flavor QCD has been successful to 
reduce the sea quark mass down to a few MeV. 
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1. Introduction 

Lattice QCD, has been successful to analyze the low-energy dynamics of 
mesons and baryons. Non-perturbative quantities, such as the hadron spec- 
trum, matrix elements, the chiral phase transition etc. have been investi- 
gated by large-scale calculations often using supercomputers. 

The lattice regularization, however, violates some important symmetries 
that the continuum theory has. For instance, the translational invariance is 
violated except for its discrete subgroup. The chiral symmetry, one of the 
most important symmetries of QCD, is also difficult to realize on the lattice. 
This is known as the fermion-doubling problem that any chiral Dirac opera- 
tor must have unphysical poles. 1 ' 2 In order to eliminate the doubler modes, 
one has to give up the chiral symmetry, or to employ complicated Dirac 
operators satisfying the Ginsparg- Wilson relation, 3-8 for which the locality 
is less obvious and more importantly much larger numerical cost is required 
compared to the explicitly local (the so-called ultra-local) Dirac operators. 
Most of the previous QCD simulations with the dynamical quarks were, 
therefore, limited to those with the Dirac operator which explicitly breaks 
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the chiral symmetry. Moreover, the sea quark mass in such simulations were 
much larger than the physical values. Their results may contain system- 
atic effects due to additive quark mass renormalizaiton, unwanted operator 
mixings with opposite chirality, additional symmetry breaking terms in the 
chiral perturbation theory, chiral extrapolation from rather heavy quark 
masses, and so on. 

The JLQCD collaboration started a new project to simulate QCD near 
the chiral limit. It uses a new supercomputer system installed at KEK in 
early 2006. We employ the overlap-Dirac operator, 6,7 which satisfies the 
Ginsparg- Wilson relation 8 and thus realizes the exact chiral symmetry. 9 In 
order to avoid gauge configurations that are too rough and the topology is 
not well-defined, we use the Iwasaki action 10 ' 11 combined with the topology 
conserving action 12-14 for the gauge part of the action. It turned out that 
keeping topology is also helpful to reduce numerical costs of the dynamical 
overlap fcrmions. On a 16 3 x 32 lattice with the lattice spacing a ~ 0.11- 
0.12 fm, we have simulated two-flavor QCD with the quark mass down to <~ 
3 MeV, which is even smaller than the physical up and down quark masses. 
We found that the locality of the overlap-Dirac operator is good enough 
compared to the QCD scale (See Ref. [15]). The preliminary results from 
the project have been reported in Refs. [16-18]. 

The outline of this article is as follows. In Sec. 2, We give a brief review 
about how the exact chiral symmetry and its topological properties can be 
realized on the lattice. The results of numerical simulations are reported 
in Sec. 3. Since our simulations are limited in a fixed topological sector 
on a fixed volume lattice ~ 2 3 x 4 fm 4 , the effects from finite volume and 
topology could be significant after reaching the chiral limit. We discuss this 
in Sec. 4. Summary and discussion are given in Sec. 5. 

2. Chiral symmetry and topology 

The difficulty of chiral symmetry on the lattice is summarized by the 
Nielsen-Ninomiya theorem: 1,2 any local Dirac operator which satisfies 
Z?75 + 75 -D = must have unphysical poles. One can easily see this in 
free fermions. By the discretization, the continuum Dirac operator in the 
momentum space, D = i^^p^, is replaced by the lattice counterpart 

^naivc = ^ + „ i7M8in(aJv) . (1) 

It has unphysical poles at p^ = n/a. Here, and d* denote the forward 
and backward subtraction, respectively. Wilson's prescription 19 to avoid 
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the unphysical poles is to add a higher derivative term 

D w = ^(d, + d;) fad; ~ 7m sin(ap„) + ^ ]T ^(p,a/2). (2) 

n 

This additional term, known as the Wilson term, gives the doublers a mass 
of the cut-off scale and let them decoupled from the theory, but it explicitly 
breaks the chiral symmetry. Since the Wilson term contains the covariant 
derivative, the eigenvalue distribution of the Dirac operator is largely de- 
formed from the continuum limit. As a result, one loses the identification 
of the zero-modes, and thus the clear definition of the quark mass. 
Neuberger's overlap-Dirac operator 6,7 is defined by 

1 / (iHw \ a 

D=-\l + 75 . =s , a = ——,aH w = j 5 (aD w - 1 - s), (3) 
a \ \/a 2 H^ J 1 + s 

where a constant s is taken in a range < s < 1. It satisfies the Ginsparg- 
Wilson relation, 8 75D + D75 = aD^^D, and the fermion action Sf — 
^2 X , 4){x)Dil){x), is exactly invariant under the chiral rotation, 9 

^e M *V, ^^e lai \ 75 = 75(1 - aD), (4) 
at finite lattice spacings. Since 1 75 ^ H 2 ™ 2 1 2 = 1, the eigenvalues of the over- 
lap Dirac operator are distributed on a circle with a radius (1 + s)/a. Note 
that this circle exactly passes the origin where we can count the number 
of the chiral zero- modes (all the zero-modes commute with 75), and the 
quark mass can be simply defined by just adding m to D. Because of the 
relation 20 

n+ - n_ = Tr7 5 /2, (5) 

where n± denotes the number of zero-modes with ± chirality, and the per- 
turbative expansion of 75, 21,22 

i 5 (x,x)/2 = ^-^tie^ pa F^ v (x)F p(r (x) + 0(a 2 ), (6) 

one can identify the index n + — n_ as the topological charge. 

One should note, however, that Eq.(3) is ill-defined when H\y has zero- 
modes. From Eq.(5), 

Q = n+ - n_ = Tr 7 " 5 /2 = Tr 75 (l - aD)/2 = Tr l5 aD/2 = ^Tr-^, (7) 

one can see that the topological charge Q changes when the eigenvalue of 
Hw crosses zero. In other words, the equation H\y = forms the topol- 
ogy boundaries in the configuration space. On the topology boundaries, 
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Hw = 0, the overlap-Dirac operator is not smooth and its locality is not 
obvious. 23 Near zero-modes of Hw also cause practical problems that the 
numerical cost of approximating 1 / \JH' W increases as the inverse of the 
lowest eigenvalue and the discontinuity of D requires a huge extra numeri- 
cal cost for the dynamical quarks. 

In order to achieve Hw 7^ at a finite lattice spacing (note that Hw 7^ 
is automatically satisfied in a — > limit), which is known as the admissibil- 
ity condition, 23,24 a transparent way is to add extra fields which produces 
a determinant, 12-14 

det H w /(H w + m 2 t ) (8) 

to the theory . It is nothing but the determinant of two-flavor Wilson 
fermions and twisted-mass ghosts with a twisted mass m t . Both of the 
additional fields have a cut-off scale mass and do not affect the low-energy 
physics. With this determinant, the overlap-Dirac operator D is smooth, 
and its locality is guaranteed. Furthermore, the numerical costs are largely 
reduced. Since this determinant does not allow the topology change during 
the hybrid Monte Carlo updates, we call it the topology conserving action. 

3. Lattice simulation 

We generate the gauge configurations of two-flavor QCD according to the 
Boltzmann weight 

exp(-/?S G ) det((l - am/2)D + mf det H W /(H W + m 2 ), (9) 

where (3Sg denotes the Iwasaki gauge action 10,11 with the coupling (5 = 2.3 
and 2.35, and (1 — am/2)D + m is an expression of the massive overlap 
fcrmion. We take a = a/1.6 and am t = 0.2 for all the simulations. On 
a 16 3 x 32 lattice, we have performed 0(1000) trajectories of the hybrid 
Monte Carlo updates in the Q = topological sector. The lattice spacing 
is estimated to be a <~ 0.11-0.12 fm from the Sommer scale r (= 0.49 fm). 
Details are seen in Refs. [16-18] 

The new supercomputer system at KEK and some algorithmic improvc- 
ments 28 ~ 30 enable us to reduce the quark mass down to ma = 0.002(^ 
3 MeV), which is even less than the physical value. In fact, we observe 
that the numerical cost, or the number of multiplication of Dw, does not 
depend too strongly on the sea quark mass m. However, we note that this 
insensitivity to m indicates that the lowest eigenvalue of D is larger than 
the quark mass, due to the finite volume effects, as discussed in the next 
section. 
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Fig. 1. Chiral extrapolation of pseudo-scalar (left) and vector meson masses (right). 

Meson masses are evaluated from the correlation functions measured 
with smeared source and local sink operators. Figure 1 shows the chiral 
extrapolation of the pion and the vector meson masses. The pion mass 
agrees well with the leading ChPT formula, to 2 cx toE / F% , where the bare 
quark mass to needs no additive renormalization, owing to the exact chiral 
symmetry. With our present statistics, we find no significant deviation from 
a simple linear fit 



TO 2 = bps TO, TOy = ay + by to 2 , 



(10) 



which gives x 2 /dof < 1.0. With this linear fit, we obtain a = 0.1312(23) fm 
using the p meson mass as an input. This is consistent with the estimate 
from rp within 10% accuracy. 



4. Finite V and fixed Q effects 

Since our simulations are done on a fixed volume lattice (~ 1.8-2 fm) 4 
in a fixed topological sector, any observable may systematic errors due to 
the finite volume and fixed topology. Although the simplest solution to 
eliminate these errors is to go to larger volumes, we can take an alternative 
way that we treat finite V and Q effects using an effective theory. 

The chiral perturbation theory (ChPT) 25 ' 26 and the chiral Random Ma- 
trix theory (ChRMT) 27 are valid in estimating the finite V effects on pions, 
especially when the quark mass is so small that the pion Compton wave 
length is larger than the lattice size. This set-up is known as the e-regime. 
In the e-regime, ChRMT describes the eigenvalue spectrum of the Dirac 
operator with the chiral condensate E as an unique parameter. The pion 
correlator in the e-regime is no more exponential but quadratic function of 
t, of which curvature depends on the pion decay constant F v . 

We compare our numerical result at the lightest quark mass with the 
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Fig. 2. Left: the cumulative distribution of the Dirac eigenvalues. Right: The pion 
correlator. 

prediction of ChRMT and ChPT. The left panel of Fig. 2 shows the cumu- 
lative eigenvalue distribution of the overlap Dirac operator, from which we 
extract the chiral condensate as £ = (251(7)(11) MeV) 3 , where the errors 
are statistical and an estimate of the higher order effects in the e-expansion. 
This value is consistent with earlier works 31,32 which are done with heavier 
quark masses. The right panel shows the pion correlator and the quadratic 
fit curve, which gives F n — 86(7) MeV. Note that this value is obtained 
near the chiral limit without doing any chiral extrapolations. 

5. Summary and discussion 

The QCD simulation in the chiral regime is feasible with the exact chiral 
symmetry respected. The topology conserving action is very helpful to re- 
duce the numerical cost of the dynamical overlap quarks. We can reduce 
the sea quark mass to the physical up and down quark masses or even 
lower. Near the chiral limit, the finite V and fixed Q effects are important 
since the pion is sensitive to these effects. Through the chiral perturbation 
theory or the chiral Random Matrix Theory, these effects can even be used 
to extract the low energy constants, such as the chiral condensate and the 
pion decay constant. To extend our lattice size is, however, important to 
confirm them in the future works. 

Numerical simulations are performed on Hitachi SRI 1000 and IBM Sys- 
tem Blue Gene Solution at High Energy Accelerator Research Organization 
(KEK) under a support of its Large Scale Simulation Program (No. 06-13). 
This work is supported in part by the Grant-in- Aid of the Ministry of Ed- 
ucation (No. 18840045). 
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